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Abstract. In this paper, two sharp inequahties for bounding the psi function 
i/i and the harmonic numbers Hn are estabhshed respectively, some results in [I. 
Muqattash and M. Yahdi, Infinite family of approximations of the Digamma 
function, Math. Comput. Modelling 43 (2006), 1329-1336.] arc improved, 
and some remarks are given. 



1. Introduction 
It ie well-known that the classical Euler's gamma function is defined by 

/"OO 

T{x) = / e-H^-^t (1) 
Jo 

for X > and the derivative of its logarithm is called the psi or digamma function 
and denoted by ip{x) for a; > 0. 

In [6], an infinite family of approximations for the psi function ipix) on (0, oo), 
denoted as {la, a e [0, 1]}, where 

Ia{x)^\n{x + a)--, (2) 

X 

was constructed. Among other things, Corollary 2.3 and Theorem 3.2 in [6] may 
be recited as follows: 

(1) For all X e (0,oo), 

\n(x + l)-->tp(x)>\nx~-. (3) 

X X 

(2) For every x e (0, oo), there exists an a G [0, 1] such that ipix) = Ia{x). 
The inequality (3) may be rearranged as 

1 > exp (^{x) + ^ - x = Q{x)>{) (4) 

for X G (0, oo). Since r(a; + 1) = xT{x) for x > Q, taking the logarithm of this 
recurrent formula and differentiating yields 

^[x + l)^'ip{x) + -. (5) 

X 

As a result, the function Q{x) defined in (4) may be rearranged as 

Q(x) = e'^("+i) - X (6) 
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for X G (0,00). The graph of Q{x) in the interval (0,9), plotted by the famous 
software Mathematica 5.2, and the limits 

lim Q{x) = and lim Q{x) = -, (7) 

2:^0+ X^QO 2 

calculated also by Mathematica 5.2, see Figure 1, show that the function Q{x) 




Figure 1. The graph on (0, 9) and two limits as a; — > 0+ and 00 of Q{x) 



is possibly decreasing on (0, 00) and \ < Q{x) < 



0.56- •• for X e (0,00), 



where 7 = 0.577 • • • stands for the Euler-Mascheroni's constant. These analyses 
imply that conclusions obtained in [6] may possibly be refined and restated more 
accurately. 

Our main results are included in the following theorems. 



Theorem 1. For all x G (0, 00), 



ln(a; + -)--< tp(x) < Inix + e '')--■ 
2 J x X 



(8) 



The constants | and e = 0.56 ■ ■■ in (8) are the best possible. 



Theorem 2. The function Q{x) defined in (6) is strictly decreasing and strictly 

convex on (— l,oo). 

Theorem 3. For every x G (0, 00), there exists a unique number a € (s"''') 5) 
such that i^{x) = Ia(x). Conversely, for every a G (e^"', 5), there exists a unique 
number x G (0, 00) such that i}{x) = Ia{x). 

In [6, Definition 3.5], the so-called error of the approximation ^'(a;) ~ Ia{x) was 
defined by 

Ea{x) = ^{x) - Ia{x) = ipix) - ln{x + a) + - (9) 

X 

for x > and a G [0, 1]. In [6, Theorem 3.7], it was proved that the errors Ea{x) 
for X e [2, 00) and a € [0, 1] are uniformly bounded between — In | and In |. This 
can be restated as follows. 

Theorem 4. Let c — 1.4616321 ■ ■ ■ is the only positive root of the psi function tp{x) 
on (0, 00). For x <E (c, 00) and a € (e~^, 5) , the following conclusions are valid: 



(1) The errors Ea{x) are uniformly bounded between — In 

(2) ^{x)=\n{x + a)-^ + 0{\n^). 



c+e 



S+T72 



and In 



c+e 



B+172' 
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It is well-known that the n-th harmonic numbers are defined for n € N by 

n 
k=l 

and that Hn can be expressed in terms of the psi function ip{x) by 

F„=V(n + l)+7- (11) 

By virtue of the decreasing monotonicity of the function Q{x) and the for- 
mula (11), the following new bounds for Hn are derived as follows. 

Theorem 5. For n e N, 



(12) 



In^n + + 7 < Hn{n) < ln(n + e^'^ - l) + 7- 
In the final section, some remarks about above conclusions are given. 

2. Proofs of theorems 

Proof of Theorem 1. Since = —7, the first limit in (7) is valid clearly. 
In [10], it was derived that 

^ <\i}{x + l)-\nx < ^ (13) 



2x 12a;2 ^ 2x 

for a; > 0. Then 

^gl/(2.)-l/(12.^) _ ^ ^ < ^^1/(2.) _ ^ 

It is easy to check that both bounds for Q{x) above tend to | as a; — > 00. The 
second limit in (7) follows. 

In [4, Lemma 1.1] and [5, Lemma 1.1], the inequality 

V''(a;)e'''(^) < 1 (14) 

for a; > was obtained. This means that 

[Q{x)]' = ij'ix + l)e'''("+i) - 1< (15) 

for x > —1. Consequently, the function Q{x) is strictly decreasing on (— l,oo). 

Combining the decreasing monotonicity with two limits in (7) of Q{x) leads 
obviously to Theorem 1. □ 

Proof of Theorem 2. The decreasing monotonicity has been proved in the proof of 
Theorem 1. 

Differentiating (15) once again gives 

[Q{x)]" = {V'"(x + 1) + [^'{x + l)]2}e'^(-+i). 
In [2, p. 208], [4, Lemma 1.1], [5, Lemma 1.1] and [7, 8], the function 

[^'{x)r+r{x) (16) 

for X € (0, 00) was verified to be positive. Hence, the function [(5(a;)]" > and 
Q{x) is strictly convex on (—1, 00). □ 

Proof of Theorem 3. Since the function Q{x) is strictly decreasing from (0, 00) onto 
, 5), that is, the mapping Q : (0, 00) 1— > 5) is bijective, then the proof of 

Theorem 3 is easily completed. □ 
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Proof of Theorem 4- By the inequaUty (8), it is easy to see that 

0<mx)-la{x)\ < |/i/2-/e-.| 

for X e (c, oo) and a £ 5), which is equivalent to 
0<|KW|<ln^^— -<ln. 



x+1/2 - c+ 1/2' 
The proof of Theorem 4 is complete. □ 

Proof of Theorem 5. Since Q{x) is strictly decreasing on (0, 00), it follows that 
lima;^oo Q{x) < Q{x) < (5(1) for x G [1, 00), which is equivalent to 

1 < eV'(-+i) -x< e*(2) _ 1 
2 

In^x + < + 1) < ln(a; + e^"'' - l). 

Taking x = n G N and using the formula (11) in the above inequality leads to the 
inequality (12). The proof of Theorem 5 is complete. □ 

3. Remarks 

Remark 1 . Note that the inequality (8) in Theorem 1 refines the double inequality 
in Corollary 2.3 of [6]. 

Remark 2. The conclusions in Theorem 3 clarify the uncertain claims between lines 
6-7 in [6, p. 1332]. 

Remark 3. In [10], the following sharp bounds for ff„ were established: For n € N, 

^'^" + ^+ 2n+l/(l-7)-2 ^^^"^<^"" + ^+^;^- ^''^ 

The constants — 2 and | are the best possible. 

In [3, pp. 386-387] and [9], alternative sharp bounds for iJ„ were presented: For 
n e N, 

1 + ln(^/i - 1) - ln(ei/("+i) - l) < if„ < 7 - ln(ei/("+i) - l). (18) 

The constants 1 + ln(Y^ — l) and 7 in (18) are the best possible. 

There have been a lot of literature devoted to bounding harmonic numbers Hn- 
For more information on Hn, please refer to [3, 9, 10] and related references therein. 

Now we compare analytically the bounds among (12), (17) and (18). For this 
purpose, let 



m 

and 



ln(^x + 1^ + 7 - [1 + ln(^/i - 1) - ln(eV(-+i) - l)] 



g{x) = \n{x + e^"^ - l) + ln(ei/(^+i) - l) . 
for a; > 1. The function f{x) may be rearranged as 



fix) = In 



+ I) (e^/(^+^) - 1)] + [7 - 1 - HVe -l)]^lnll^+ A, 
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where 



(2a;2 + 2X + 1) [eV(^+i) - 2{x + l)^/{2x^ + 2x + l) 



[x + 1)2 

^ (2a;2 + 2a;+l)/2(l/(a; + l)) 



f2{u) 



[l + {l-uf 



1 + (1-U)2 1 + l + (l_u)2 

and /3(u) = M^e" > for u e (0,5]. Since /3('u) is increasing on (O, ^] and 

lim„^o+ fsi^) = 0, it follows that f3{u) > and /2(m) > on (O, ^] . This leads 
to f{{x) > and fi{x) being strictly increasing for a; > 1. Therefore, the function 
f{x) is increasing for a; > 1, with 

/(I) =ln[(^/i - l)/2] + 7- 1 = -0.01731922699030 •• • 

and 

+ 1/2 ei/(^+i)-l" 



lim /(x) 



lim In 

x—>-oo 



+ 7- 1 -ln(Ve - 1) 



x+1 l/{x + l) 
= 7- l-ln(\/i - 1) 
= 0.00996779446872 • • • . 

As a result, the left-hand side inequalities (12) and (18) are not included each other. 
Numerical calculation gives 

/(2) = -0.001061745178- •• and /(3) = 0.004039213518 ••• , 

thus, as n > 3, the lower bound in (12) is better than the one in (18). 
Numerical computation gives 

g{2) = -0.00060205073286 • • • and g{3) = 0.00153070402207 • • • . 

It is easy to see that 

"a; + 6^-^-1 6^/(^+1) -1" 



lim g{x) = In 



0. 



x + 1 l/(a;+l) 
Combining with the graph of g{x) on (1, 29), see Figure 2, it can be possibly claimed 
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Figure 2. The graph of g{x) on the interval (1, 29) 

that, as n > 3, the upper bound in (12) is not better than the one in (18). 

Similarly, it may be showed that the upper bound in (12) is not better than the 
one in (17) as n > 2, that the lower bound in (12) is not better than the one in (17) 
at all, and that the inequality (18) is not better than the one (17) at all. Thus, 
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it may be said that the inequahty (17) is the best among inequahties (12), (17) 
and (18). 

Remark 4. In [7, 8], the positivity of the function (16) was generaUzed to a much 
general result including, as a particular case, its being completely monotonic on 

(0,c^). 

Remark 5. It is conjectured that the function Q{x) defined in (6) is completely 
monotonic on (— l,oo), that is, {—l)'^[Q{x)]^'''> > on (— l,oo) for A; > 0. 
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